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Resumen 


La composición anual de la evaporación en la cuenca del lago de Pátzcuaro se 
integró por la evapotranspiración de la vegetación, las pérdidas por intercepción, 
la evaporación del cuerpo de agua y la evapotranspiración de la vegetación 
hidrófita. La evapotranspiración de la vegetación se estimó con los registros de seis 
estaciones climatológicas clásicas y siete evaporímetros Elgage. La evaporación 
del lago se midió con el método BREB. La evaporación por la intercepción de la 
vegetación fue calculada con relaciones lineales en función de la precipitación para 
diversos tipos de vegetación y la transpiración de la vegetación hidrófita a partir de 
mediciones puntuales con una cámara de circuito abierto. La evapotranspiración 
resultó de 308 mm < ETR < 503 mm (+ 60), la evaporación del lago de 1 179 
mm, la evapotranspiración de la vegetación hidrófita (ETH) de 1 824 mm y la 
correspondiente a la intercepción entre 102 mm < El < 227 mm (+ 51). En el contexto 
del balance hídrico, en el nivel subcuenca se encontró un déficit de 12 mm en la de 
Ajuno y más del 65% en la propia cuenca del lago. En el resto de las subcuencas, el 
balance del residuo fue positivo para la generación de escurrimientos superficiales 
e infiltración, y/o percolación profunda, entre 18 mm < (P-ET) < 156 mm (+ 53). Los 
resultados encontrados muestran la importancia de la vegetación y el lago en la 
distribución de la evaporación de la cuenca. 
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Introducción 


La evaporación generada por los elementos 
que constituyen la superficie de una cuenca, 
la transpiración de la cobertura vegetal, la 
evaporación del suelo desnudo, la evaporación 
de los cuerpos de agua, la evaporación generada 
por la intercepción de lluvia y la producida por 
la transpiración de la vegetación hidrófita son 
componentes primordiales del balance hídrico. 
Los cambios en la proporción de la evaporación 
a largo plazo pueden relacionarse con las 


modificaciones en el uso del suelo, la cobertura 
vegetal y el clima (Bouwer et al., 2008). 

para 
la evaporación se han originado en climas 


Numerosos métodos determinar 
contrastantes y diferentes principios teóricos, 
los cuales pocas veces son comprobados con 
mediciones directas de la evaporación cuando 
se aplican a otros climas o condiciones 
diferentes a las de su origen. El principal 
problema se debe a la disponibilidad y calidad 
de observaciones directas de la evaporación, 
sobre todo para los cuerpos de agua. 
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En estudios de procesos hidrológicos y 
ecológicos, 
aplicación del método está sujeta al tipo de 


inevitablemente la elección y 
información disponible. Son ampliamente 
utilizados por su simplicidad los métodos: 
1) Pan (tanque clase A) y los clasificados 
dentro del grupo de temperatura: Papadakis 
y Thornthwaite; 2) los que combinan la 
temperatura y duración del día: Blaney-Criddle 
y Hamon (Rosenberry et al., 2007); 3) los que 
requieren información más detallada y cuentan 
con bases teóricas consistentes son clasificados 
dentro del grupo de Dalton: transferencia de 
masa (Harbeck et al., 1958) y Ryan-Herleman; 
4) o bien los que consideran la radiación solar 
y la temperatura: Makkink, Jensen-Haise y 
Sthepens-Stewart; 5) los que necesitan más 
información son llamados métodos combina- 
deBruin, 

Penman, Brutsaert-Stricker y deBruin. 


dos: —Priestley-Taylor, Keijman, 

En contraste, el método BREB (Bowen 
Ratio Energy Balance) es un procedimiento 
consistente e incluso un estándar en Estados 
Unidos por su confiabilidad y robustez; sin 
embargo, por su formulación, el método BREB 
ignora cualquier covarianza entre la velocidad 
del viento y los gradientes de temperatura 
y presión de vapor, lo cual podría provocar 
errores adicionales en el método. 

En cambio, el procedimiento Eddy Cova- 
riance (EC), a pesar de ser una técnica más 
directa para medir la evaporación, requiere 
de instrumentos costosos, limitando su uso 
en forma generalizada. Si bien algunos de los 
métodos indicados son aplicables a lagos de 
montaña (Ronsenberry et al., 2007), la decisión 
en la selección del método depende igualmente 
de la escala del tiempo de estimación de la 
evaporación, del tiempo de registro de las 
variables disponibles; es decir, su validez está 
sujeta a las condiciones particulares de cada 
sitio. 

De la misma forma, el conocimiento de la 
evaporación generada en los lagos, humedales 
O bien presas es un componente valioso en 
la estimación del balance y la gestión de los 
recursos hídricos, considerando sus diversos 


usos dentro del sector urbano, industrial y de 
producción de energía eléctrica. 

La certeza en la estimación de la evapora- 
ción de los lagos asegura en consecuencia 
una buena planeación en el ahorro del agua 
de los sistemas de abastecimiento. Por tal 
razón, el disponer de mediciones directas de 
la evaporación de la cuenca del lago de Pátz- 
cuaro es parte fundamental para un adecuado 
hídrico de la 
conservando en todo momento los aspectos 


aprovechamiento región, 
ecológicos a corto y largo plazos. 

En el país, la Comisión Nacional del Agua 
(Conagua) reportó para el periodo 1941-2000 
una precipitación normal en el territorio de 
771 mm, donde el 67.3% se concentra en cuatro 
meses (junio-septiembre). Específicamente 
en el estado de Michoacán, la precipitación 
fue de 802 mm. Maderey y Torres-Ruata 
(1990) estimaron una evaporación entre 800 
y 900 mm por año en la región del lago de 
Pátzcuaro, aplicando el método de Turc (1954), 
basándose en la información de 543 estaciones 
climatológicas, con registros del periodo 1945- 
1980. 

El residuo de la precipitación y la 
evaporación revela que los escurrimientos 
superficiales einfiltración y / orecarga profunda 
son prácticamente nulos. Cabe señalar que 
no obstante que el método de Turc (1954) es 
ampliamente utilizado por su simplicidad, no 
estima adecuadamente la evapotranspiración 
anual. Incluso a escala decenal, Seguin (1975) 
demuestra que subestima en un 12% la 
evapotranspiración real. 

El Instituto Mexicano de Tecnología 
del Agua (IMTA) llevó a cabo en 2004 un 
estudio del balance hídrico en la cuenca del 
lago de Pátzcuaro (Lafragua y Gutiérrez, 
2005) integrándolo en tres flujos: superficial, 
subterráneo y del lago. El balance hídrico 
mostró que el 77% de la precipitación 
regresa a la atmósfera en forma de vapor, el 
10% se infiltra y únicamente el 13% genera 
escurrimientos superficiales. Los flujos fueron 
estimados con la base de datos del año 2004, 
generada con una red de cinco pluviómetros 


González-Sosa, E. et al., La evaporación en la cuenca del lago de Pátzcuaro, México 


digitales y tres estaciones meteorológicas 
clásicas automatizadas (EMASs). 

Las EMAs 
variables climatológicas y la precipitación cada 


estaciones registraron las 
diez minutos; en cambio, los pluviómetros 
La 


evaporación total de la cuenca incluyó la 


registraron cada treinta minutos. 
producida por el lago, la vegetación de la 
cuenca, la vegetación hidrófita y las pérdidas 
por intercepción. La combinación del modelo 
de Zhang et al. (2001) y el criterio de Turc 
(1954), en su modalidad anual, produjeron 
por parte del sistema suelo-planta una 
evapotranspiración de 1 193.85 mm, superior 
en 1.28 veces a la precipitación media anual 
(931.88 mm); es decir, un cociente de ETP/P 
por arriba del valor medio de 0.65 referido 
por diversos autores (e.g. Brutsaert, 1982, 
1991). Por subcuenca, dicha relación osciló 1.07 
< ETP/P < 1.45, donde los valores extremos 
ocurrieron al sur y norte del lago, Pátzcuaro y 
San Jerónimo, respectivamente. Las pérdidas 
por 
como el 18% de la precipitación anual, 


intercepción (El) fueron estimadas 
admitiendo implícitamente una capacidad 
de almacenamiento y distribución uniforme 
de la vegetación en 41.3% [bosque-(330.74 
km2)] del 80% de la cuenca. De esta forma, la 
integración de las superficies boscosas de las 
subcuencas del sistema del lago de Pátzcuaro 
generaron una lámina de El de 174.76 mm 
(18.7%). En tanto, con la base de datos de la 
estación La Pacanda —situada en la parte alta 
de la isla del mismo nombre, 30 m por arriba 
de la superficie del agua— se estimó el flujo 
de vapor del lago por separado de los otros 
flujos evaporativos de la cuenca a través de 
la ecuación de Penman (1948), sin considerar 
las superficies de las islas y la vegetación 
hidrófita. El flujo de vapor resultó de 989.5 
mm, superior en 1.12 veces a la lámina de 
lluvia (879 mm) sobre la superficie del lago. 
La transpiración de la vegetación hidrófita 
fue deducida en proporción lineal al flujo 
de vapor del cuerpo de agua del lago, como 
1.25 veces la evaporación anual, equivalente a 
274.25 mm (Lafragua y Gutiérrez, 2005). 


El objetivo de este trabajo fue evaluar la 
evaporación en la cuenca del lago de Pátzcuaro 
a partir de mediciones directas y patentizar 
la importancia de la evaporación del lago 
y la producida por la vegetación terrestre e 
hidrófita, así como establecer una comparación 
entre la evapotranspiración calculada con el 
método de Penman-Monteith-FAO (Allen 
et al., 1998) y la observada con el ETgage, 
comparar la evaporación del lago estimada 
con la ecuación de Penman (1948) y la medida 
con el método BREB, y conocer la distribución 
de la evaporación generada por los diversos 
componentes de la cuenca para establecer el 
balance hídrico. 


Metodología 


De acuerdo con los criterios propuestos por 
Ahuja (1959), Karasseft (1986) y Schaake 
(2000), correspondientes a los criterios de 
la Organización Meteorológica Mundial 
(OMM), se evaluó la densidad de las estaciones 
para la cuenca del lago de Pátzcuaro. Con base 
en el estudio realizado por Lafragua y Gutiérrez 
(2005), fueron seleccionados siete sitios para 
medir la evapotranspiración con atmómetros 
(ETgage modelo E, cubierta 4G2, que impide 
el paso de la precipitación y permite la 
evaporación durante la lluvia; Berrada et al., 
2001), con el objetivo de cubrir siete subcuencas 
de la cuenca. Los ETgage fueron instalados a 1.0 
m de altura en áreas libres de obstáculos para 
evitar errores por el paso del sol y del viento 
(Manual ETgage). 

Las figuras 1 y 2 muestran la ubicación 
en planta y en altitud de las estaciones y 
ETgage. Para evaluar la ET'P con la ecuación 
de Penman-Monteith-FAO (Allen et al., 1998) 
se instalaron tres estaciones meteorológicas 
automáticas tipo WATCHDOG-900, adicio- 
nales a las tres estaciones EMAs (IMTA), 
cumpliendo con ello la recomendación 
internacional para la densidad de redes 
de estaciones climatológicas (OMM, 1970, 
1994). Una en la subcuenca de Ajuno, otra 
en la subcuenca de Tzurumútaro y la tercera 
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en la zona de San Andrés. Los ETgage y las 
estaciones de los sitios de Pichátaro y Yotatiro, 
en la zona poniente de la cuenca a 2 600 y 2 
200 msmm, respectivamente. 

Para medir la evaporación del lago se ins- 
taló una estación climática Bowen (Campbell 
Scientific, 1998) basada en el método de 
BREB, que mide los componentes del balance 
de energía de superficie y los gradientes de 
temperatura y presión de vapor del aire. La 
estación fue colocada en la vertiente norte de 
la parte central del lago, al pie de la isla de 
La Pacanda, con el propósito de garantizar un 
fetch conveniente (F > 150 m). Los gradientes 
de presión de vapor y de temperatura se 
midieron con una separación de 1 m, con 
un desnivel de 3 m respecto a la superficie 
del agua y registro cada veinte minutos. La 
temperatura del aire (T), la humedad relativa 
(H), y la velocidad (U,) y dirección del viento 
fueron medidas a una altura 4 m por arriba del 
nivel del agua. Durante dos meses se midió 
la temperatura del agua con un termopar de 
cobre suspendido a 2 cm de profundidad. 

Adicionalmente se llevaron a cabo dos 
campañas de medición de la transpiración 
de la vegetación hidrófita mediante una 
cámara portátil de circuito abierto, la cual se 
basa en la medición del flujo de vapor que 
circula al interior de la cámara (Stannard y 
Weltz, 2006). A través de imágenes de satélite 
y verificación sobre el terreno se delimitó 
la superficie y composición de las especies 
dominantes de la vegetación tanto terrestre 
como hidrófita. 


Método Bowen Ratio Energy Budget (BREB) 


La relación de Bowen (PB) representa el cociente 
entre los flujos de calor latente (H) y el flujo 
de calor sensible (LE). La relación de Bowen, 
combinada con el balance de energía de 
superficie, definen la formulación del método 
Bowen Ratio Energy Budget (BREB) (Stannard, 
1997): 


H 
E (1) 


R, =H+LE+G (2) 


donde R, es la radiación neta y L es el calor 
latente de vaporización del agua (2450 J g'). 
La combinación de las ecuaciones (1) y (2) 
producen la ecuación para estimar el flujo de 
vapor: 

1 


LE==— 
1+p 


(R, -G) (3) 
Dentro de un periodo promedio (30 min) 

se han formulado relaciones empíricas que 

relacionan los flujos y gradientes verticales: 


¿€ 
pe pj de 
Y Oz 
se (4) 
H=-p¿Cyky — 
Pa p*h Oz 


donde C, es el calor volumétrico del aire; e, la 
presión parcial de vapor; p,, la densidad del 
aire; k, el coeficiente de difusión de vapor, y k, 
esel coeficiente de difusión de calor. Asumiendo 
k, = k, y disponiendo de los gradientes en 


dos niveles dentro de una superficie dada, la 
relación de Bowen resulta: 


oT 

_ 9 _, AT 

Pe Y de (5) 
Oz 


Siendo y la constante psicométrica. Para el 
caso de cuerpos de agua, la ecuación (3) puede 
definirse de la forma siguiente (Rosenberry et 
al., 2007): 


1 
= == = 6 
E olLa+prrcr] Q:+Q.-Qu) (6) 


siendo C el calor específico del agua [4 186 
J kg* *C"]; R,, la radiación neta [Wm>]; T,, 
la temperatura de la superficie del agua 
[*C]; Q, el cambio de calor almacenado en 
el cuerpo de agua [Wm-”]; Q,, la energía 
neta (positiva cuando es por advección del 
algo) a partir de la precipitación, superficie 
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del agua y acuífero [Wm”], y Q, la energía 
neta por conducción de los sedimentos en el 
lago [Wm-”]. La relación de Bowen toma la 
forma: 


E=T 
PB = CP Per (7) 


eo a 


C, es el coeficiente empírico determinado 
por Bowen igual 0.61 (*C*); P, la presión 
atmosférica (kPa); T, la temperatura del aire 
(*C); T, la temperatura de la superficie del 
agua; €, la presión de vapor a saturación para 
la temperatura de la superficie del agua (Pa), y 
e, es la presión de vapor (Pa). 

Como la profundidad actual del lago es 
somera, entre 2 y 12.5 m, y no fue medido el 
perfil de temperatura, se estableció la hipótesis 
de una sola termoclina. Al igual que el trabajo 
de Rosenberry et al. (2007), quienes reportan 
que los valores de Q, y Q, generalmente son 
pequeños, en este estudio igualmente son 
ignorados tales valores. 

Para evaluar la bondad del método de BREB 
se comparó la evaporación observada con la 
estimada con el método de Penman (1948). Sin 
duda por la dificultad natural de disponer las 
mediciones convenientes para la evaluación 
de la evaporación del lago de Pátzcuaro fue 
un desafío establecer el sitio experimental 
y garantizar un buen fetch, máxime que las 
mediciones deben ser representativas del flujo 
de vapor generado por el lago. 

El sector norte presentó una distancia 
suficientemente larga como para suponer 
que las mediciones están dentro de la 
subcapa atmosférica de equilibrio, porción 
de la atmósfera en equilibrio con la superficie 
evaporante. Como los requerimientos del fetch 
también están asociados con la rugosidad y 
condiciones atmosféricas del sitio, se siguió la 
regla 1:100 de la relación altura-longitud para 
garantizar un grado razonable de equilibrio 
entre el sensor y la superficie evaporante 
(Heilman ef al., 1989). 


Descripción del sitio 


El lago de Pátzcuaro se localiza dentro 
de la cuenca del mismo nombre, entre las 
coordenadas 19? 27” a 19” 44” de latitud norte 
y entre los 101? 26' a 101? 53” de longitud 
oeste, parte central del estado de Michoacán, 
cuenca Lerma-Santiago. La cuenca es del tipo 
endorreica y se sitúa dentro de lasubprovincia 
neovolcánica tarasca perteneciente al eje 
Neovolcánico Transversal. El clima en la 
cuenca es templado subhúmedo, con lluvias 
de verano. Una temperatura media anual 
14 *C <T<20*C y una precipitación media 
anual histórica del periodo 1963-2003 entre 
900 mm < P<1 140 mm. La geomorfología de 
la cuenca da origen a cinco climas, 18 tipos de 
suelos y tres tipos de bosques mixtos: pino- 
encino, de pino (Pinus), encino (Quercus) y de 
oyameles (Abies religiosa), y cuatro especies 
arbustivas, incluso asociaciones de huizache- 
mezquite. 

La agricultura ocupa el 40% de la super- 
ficie de la cuenca y los cultivos se ubican en 
la ribera del lago y la zona de piedemonte. En 
la figura 1 se muestran los tipos dominantes 
de vegetación obtenidos de una imagen 
de satélite del año 2007. Puede observarse 
que en los bordes del lago se tiene la 
formación de vegetación hidrófita, acuática 
y subacuática, dominada por tules (Scirpus 
americanus, Typha latifolia, Sagitaria gramínea, 
Cyperus niger). 

El lago de Pátzcuaro puede considerarse 
como un lago tropical de montaña, con una 
profundidad somera y variable debido a su 
topografía y la gran dinámica del proceso 
erosión-sedimentación; la longitud y anchura 
máxima son de 20 y 10 km, respectivamente. 
Las profundidades varían entre 2 y 4 m 
en la fracción sur y sureste (Ihuatzio y 
Erongarícuaro), y 12.5 m en el sector oriente 
(Quiroga). En la zona denominada cuello del 
lago, donde se localizan las islas de Janitzio, 
Tecuén, Yunuén y La Pacanda, la profundidad 
es de 5 m. 
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Figura 1. a) Ubicación de los sitios de medición de la evaporación y las variables climatológicas. La ubicación 
de los sitios respecto a la altitud sobre el trazo de cuatro ejes transversales: dos en la dirección norte-sur 
(2-2” y 4-47) y dos en la dirección este-oeste (1-1” y 3-37); b) distribución de la cobertura vegetal, clasificada 
con base en la imagen de satélite de 2007. 


Ciencias del Agua, julio-septiembre de 2010 


( vj30]0u32], 


González-Sosa, E. et al., La evaporación en la cuenca del lago de Pátzcuaro, México 


Resultados y discusión 


Las variables registradas en las estaciones 
EMA's mostraron un comportamiento pare- 
cido, salvo la temperatura en Pichátaro, 
que fue inferior al resto de las estaciones, 
originado por su ubicación en la montaña alta 
de la cuenca (figuras 1 y 2). 

En la figura 3 se presentan los valores 
medios diarios de las estaciones de: Tupátaro, 
Pichátaro, Yotatiro, Chupícuaro, La Pacanda 
y Pátzcuaro, durante el periodo de julio 2006 
a junio 2007: temperatura 5% C< T,<25*C (+ 
2.1, +2.6, + 2.4, + 2.2, + 244, + 2.8); humedad 
relativa 21% < H,< 94% (+ 16.9, + 14.0, + 13.7, 
+ 12.9, + 11.8, + 13.7); déficit de presión de 
vapor 0.07 KPa < DP < 1.7 KPa (2 0.3, + 0.2, 
+ 0.3, + 0.3, + 0.3, + 0.3); velocidad del viento 
0.01 ms*<U <5.7ms* e 0,3, « 04d, 50,2, + 
0.4, + 0.7, + 0.3); la radiación global 7 Wm?< 
R,< 371 Wm?” (+ 55.6, + 56.4, + 51.3, + 56.7, + 
56.1), y una radiación neta 28 Wm”< R, < 253 
Win*-(+30.8, e 31:6, 432.1, +38:.7, 37.9) 

El comportamiento 


semejante de las 
variables climáticas en los seis sitios (figura 
3) evidencia un clima regional dentro de la 
cuenca. En contraste, los registros de la estación 
Bowen y de la isla de La Pacanda mostraron 
un comportamiento diferente al observado en 
otros sitios, un clima particular del lago. Para 
el primer semestre de 2007: temperatura 12 *C 
< T,<20*C (+ 1.7); humedad relativa 35% < H, 
< 80% (+ 58.8); déficit de presión de vapor 0.4 
KPa < DPV < 1.4 KPa (+ 0.82); velocidad del 
viento 0.5 ms*< U,< 9 ms” (+ 13.6); radiación 
global 95 Wm?< R, < 350 Wm? (+ 47.5), y 
una radiación neta 55 Wm”< R,< 250 Wm? 
(+ 41), siendo las variables la, DPV y Hr las 
que mostraron un comportamiento diferente 
respecto a las otras estaciones (figura 3). 

Los vientos dominantes circularon en 
dirección del sector noroeste de la cuenca, con 
mayor potencia en la zona central del lago, 
incluso con ráfagas instantáneas máximas del 
orden de 100 km h”. 

Tal comportamiento del viento revela 
que es un elemento regulador del clima de 


la cuenca debido a su acción convectiva- 
advectiva como fuente generadora de flujo de 
vapor y humedad en el interior de la cuenca, 
primordialmente en el sector norte del lago. 
Por otra parte, las velocidades reducidas de 
Yotatiro y Pátzcuaro, en comparación con las 
otras estaciones, hacen suponer la presencia 
de obstrucciones alrededor de los sitios, 
lo cual sin duda afecta la estimación de la 
evaporación cuando son aplicados métodos 
que consideran el viento. 

La precipitación histórica de la estación 
Pátzcuaro, periodo 1969-1991, comparada 
con la correspondiente al periodo de estudio, 
muestra que fue un año relativamente seco 
(figura 4). En el cuadro 1 se presenta un 
resumen de la distribución la precipitación 
por mes y subcuenca, así como los valores 
anuales. En el periodo de estudio, la 
precipitación anual por subcuenca fluctuó 
entre 400 mm < P < 810 mm, siendo la 
subcuenca de Ajuno con menor precipitación 
(309 mm) y la de Tzurumútaro la más 
lluviosa (809 mm), con una dispersión de 
+125 mm entre las subcuencas, incluyendo la 
propia del lago. El sector norte de la cuenca 
fue el más lluvioso, en tanto en la parte 
suroeste, zona montañosa, la precipitación 
se redujo en un 50%, y en el resto de las 
subcuencas la diferencia fue entre el 3 y 
14%, siendo significativa la diferencia entre 
el sector norte y el resto de los sectores. Las 
diferencias entre las subcuencas de Ajuno 
y Tzurumútaro pueden explicarse por la 
orografía de la zona, es decir, por la presencia 
de agua atmosférica transportada por niebla 
advectiva (Labajo et al., 1995), condensación 
típica de zonas de montaña, la cual no fue 
registrada con los pluviómetros utilizados 
en el estudio. 

También la niebla representa una signi- 
ficante fuente de vapor por la intercepción 
por condensación, la cual no es considerada 
en los métodos aplicados en este trabajo. La 
mayor parte de la precipitación se presentó 
en el ciclo verano-primavera, correspondiente 
a los meses de julio-septiembre 2006 y abril- 
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junio 2007, y sólo el 10% en el periodo otoño- 
invierno. 

La precipitación se concentró en tres 
sectores. Por orden de magnitud, el sector 
oriente, zona de piedemonte-valle, integrada 
por las subcuencas de Pátzcuaro, Tzurumútaro 
y Quiroga, resultó ser el más lluvioso, con 
810 mm año”. El sector central, integrado por 
las subcuencas de Cucuchucho, 
Tzintzuntzan, Ojo de Agua y El Lago, 


Granada, 


conjuntamente con Quiroga, conprecipitaciones 
intermedias, 730-770 mm, originadas por 
los efectos de las corrientes de humedad del 
Pacífico y del golfo. El de menor precipitación 
fue el sector norponiente, correspondiente a la 
zona montañosa, donde se ubica Ajuno, entre 
400-550 mm. En la figura 4b se presenta, a 
título de ejemplo, la distribución espacial de la 
precipitación para el mes de junio 2007. 

Cabe mencionar que análisis estadísticos 
de rupturas de series de precipitación 
muestran que existen dos grandes periodos 
históricos de lluvias. El primero de 1949 a 
1980 y el segundo de 1981 a 2004, lo que hace 
suponer que el comportamiento hasta el año 
2006-2007 continúa dentro del mismo régimen 
(Gutierrez et al., 2004). La comparación de 
la evapotranspiración observada con los 
seis Elgage situados en la zona de montaña, 
piedemonte-valle y valle, calculada con los 
métodos de Turc, Turc modificado (Abtew, 
1996, 2001), Pryestley-Taylor e Imak, Penman- 
Monteith-FAO, 
fue el adecuado para su aplicación en el 


confirmó que este último 


presente trabajo. En la figura 5a se muestra 
la evaporación mensual calculada con los 
métodos mencionados en la estación Yotatiro, 
localizada en la cuenca de Ajuno, y en la figura 
5b se tiene la evaporación correspondiente para 
el lago. 

La comparación entre el Elgage y 
los métodos también puso en evidencia 
las limitaciones de su aplicación a las 
condiciones locales de cada sitio. A título 
de ejemplo, el sitio de Pichátaro, por sus 
condiciones climatológicas, zona de lluvia 
orográfica y niebla, zona de montaña y la 


evaporación calculada con los métodos fue 
inferior a la observada con el ETgage, e 
incluso con el método de Penman-Monteith, 
a pesar de ser un método estandarizado por 
la FAO y por diversas investigaciones (e.g. 
Brown, 1981). 

En la figura 6a se ejemplifica la distribución 
espacial de la evapotranspiración calculada con 
Penman-Monteith para el mes de junio 2007. 
En la figura 6b se presenta la distribución de 
la evapotranspiración real, una vez aplicada 
la corrección de Zeng (Sharif y Miller, 2006), 
también para el mismo mes. 

La distribución espacial de la evaporación 
ocasionada por las pérdidas por intercepción 
para los diferentes tipos de cobertura no se 
muestra aquí. Al respecto, es importante desta- 
car que no son considerados exhaustivamente 
todos los tipos de vegetación dentro de las siete 
clases consideradas en el estudio. En el cuadro 
2 se resumen la distribución temporal de la 
evapotranspiración y la evaporación. 

En el cuadro 2 se muestran los valores 
mensuales de ETR corregidos con el método de 
Zeng (Sharif y Miller, 2006). Como puede verse 
a nivel potencial, los valores anuales fluctúan 
de 1 083 mm año” < ETP < 1 385 mm año”, con 
una evaporación media para la cuenca de 1 261 
(+ 97) mm año”, los cuales superan los 900 mm 
obtenidos por Maderey y Torres Ruata (1990), 
quienes aplicaron el método de Turc (1954). 
Las mayores tasas ocurren en el periodo prima- 
vera-verano, 128 mm mes” (+ 11.0) < ETP < 
170 mm mes” (+ 18.0), equivalentes a 4.3 y 5.7 
mm día”, siendo la subcuenca Ojo de Agua, 
situada en el costado interior del lago (sector 
centro) y la de Erongarícuaro (sector oeste) la 
de mayor y menor ETP, respectivamente. Una 
reducción sustancial de la ETP aparece en la 
época septiembre-enero, de 86 (+ 6.0) mm 
mes” < ETP < 119 mm mes” (+ 8), perteneciente 
a una temporada de bajas temperaturas y 
menor radiación. La aplicación del índice 
de humedad definido por Zeng et al. (2000) 
produce una ETR en un orden de magnitud 
inferior del 50% respecto a la ETP. Los valores 
anuales fluctúan de 308 mm año*< ETR < 503 
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mm año”, con una media para la cuenca de 
440 (+ 60) mm año”. En marzo-abril 2007, la 
evapotranspiración prácticamente es nula, lo 
cual coincide con un periodo de poca humedad 
relativa y un significante DPV, tal como puede 
observarse en la figura 3. Inclusive el periodo 
de sequía, que abarca seis meses, noviembre 
2000-abril 2007, con una ETR < 1.0 mm día” La 
distribución espacial ubica las mayores tasas 
de ETR en el sector sur-oeste, la cual coincide 
con la mayor zona de cultivos, con 1.2 mm día” 
< ETR < 3.4 mm día” El análisis de los errores 
de la relación de Bowen en condiciones no 
advectivas como establecen Pérez et al. (1999) 
y XiaoFeng Guo et al. (2007), entre otros, para 
el caso de superficies vegetadas, no fueron 
significativos, porque la relación de Bowen 
siempre fue mayor que -1. Cabe destacar que 


el 64% de los datos cumplieron con la relación 


Ae +pAT z 
Fha-G > (Pérez et al., 1999). 


En este mismo sentido debe señalarse 
que el transporte de flujo de vapor por 
advección en cuerpos de agua está regido 
por 
poco descritos en la literatura especializada. 


mecanismos turbulentos complejos 
En la figura 7a se presenta la comparación 
de la evaporación registrada con el ETgage 
situado en la isla de La Pacanda y el método 
de Penman (1948), a partir de la información 
de la estación climática situada también en 
la isla y la estimada con el método BREB 
(ecuación (3)). Puede observarse que existe 
una muy buena correspondencia entre el 
método BREB y el ETgage, un R?> 0.9, como 
también lo observó Puryear (2005). En tanto 
el método de Penman (1948) muestra una 
sistemática bien definida 
respecto al BREB y el ETgage (figura 7b). 
Con base en esta comparación y correlación 


subestimación 


entre Penman (1948) y Bowen se calculó la 
evaporación anual del lago. La evaporación 
del lago fluctuó 2.2 mm día*< Ew < 5.0 
mm día” y durante el periodo seco, el lago 
proporciona la humedad a la región central 
de la cuenca para satisfacer la demanda 
evaporativa de la atmósfera. Del mismo 


puede notarse un comportamiento semejante 
entre la ETR y la Ew durante el periodo de 
lluvias julio-septiembre, cuando el suelo 
alcanza la saturación y existe suficiente agua 
en la zona de las raíces. En la figura 7b se 
tienen los valores de la relación Bowen y los 
gradientes de presión de vapor y temperatura. 
La transpiración mensual de la vegetación 
hidrófita del lago se estimó con base en la 
medición en las tres especies dominantes 
(thypha, tule y lirio) y el agua, y la 
correlación con la evaporación con el método 
BREB. El cuadro 3 muestra ambos resultados: 
evaporación del lago y transpiración de la 
vegetación hidrófita. 

Puede apreciarse que la transpiración de 
la vegetación hidrófita (ETH) es persistente 
mensualmente, en comparación con la del 
lago. En el periodo agosto (2006) a enero 
(2007), la transpiración generada por la 
typha y tule superan la producida por el 
lago. Es decir, la evaporación generada por la 
vegetación hidrófita por unidad de superficie 
es significativamente mayor a la del lago, en 
un orden de magnitud de 1.3 <(ETH/Ew) <3.8. 
La evaporación producida por la intercepción 
de la vegetación, sin incluir la parte de las 
herbáceas y los lechos de residuos vegetales 
alcanza valores entre 102 mm < El < 277 mm 
(+ 51), es decir, representa por subcuenca del 
orden del 35y 19% de la precipitación anual. 
Las pérdidas por intercepción se concentran en 
el sector oeste-central del lago, Tzintzuntzan 
y Ojo de Agua, subcuencas integradas por 
cultivos y vegetación herbácea. En promedio, 
la tasa de El no rebasa 1.7 mm día”, y en el 
periodo noviembre (2006) a mayo (2007) son 
prácticamente nulas por la escasa precipitación 
en este periodo. 


Conclusiones 


En este trabajo se determinó la composición 
anual de la evaporación integrada por la 
evapotranspiración de la vegetación, las 
pérdidas por intercepción, la evaporación del 
cuerpo de agua y la evapotranspiración de la 
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Cuadro 1. Precipitación mensual en mm, julio 2006-junio 2007. 


2006 2007 
Cuenca Área | Julio | Agosto | Septiembre | Octubre | Noviembre | Diciembre | Enero | Febrero | Marzo | Abrill | Mayo | Junio | Total 
Ajuno 243 | 58 74 52 4 13 22 0 p) 19 87 | 398 
Pátzcuaro 57 | 159 | 178 114 5 17 17 0 0 38 122 | 798 
Tzurumútaro | 176 | 154 | 164 170 4 il 16 8 il 43 110 | 810 
Quiroga 125 | 141 172 213 4 8 14 3 2 16 104 | 769 
San Jerónimo cn | 165 [| 165 200 3 Y 22 0 0 7 70 | 745 
Erongarícuaro | 115 77 91 80 4 12 26 0 0 9 102 | 566 
Cucuchucho 120 [11420 157 154 4 y 17 0 0 17 84 | 697 
Granada 12 | 150 | 159 190 104 21 3 8 24 0 0 9 60 | 728 
Tzintzuntzan 120 [11580 P2U70 211 93 14 4 8 19 5) 0 16 87 | 782 
Ojo de Agua 4 | 145 | 175 216 82 11 5 7 15 3 2 13 102 | 775 
Lago 96 | 147 | 162 181 101 16 4 8 20 0 0 11 80 | 730 
Máximo 159 178 216 147 21 5 17 26 8 5 42 122 | 810 
Mínimo 58 74 52 56 8 3 Y) 14 0 0 7 60 | 398 
Media 134 | 151 160 100 14 4 11 19 Z 1 19 93 | 707 
Desv. est. 34 35 56 25 4 1 4 4 3 2 112 18 | 123 


Cuadro 2. Repartición mensual de la evapotranspiración. 


2006 2007 

Cuenca Julio | Agosto|Septiembre | Octubre | Noviembre | Diciembre | Enero | Febrero | Marzo | Abril | Mayo | Junio Total 

Ajuno 107/46| 100/54| 96/41 94/43 74/8 81/4 83/12 | 76/19 | 106/0|119/4| 116/18 | 96/59 | 1148/308 
Pátzcuaro 112/92| 101/95| 98/72 90/77 80/14 78/5 78/16| 96/16 | 136/0 |152/0| 145/33 |128/83 | 1294/503 
Tzurumútaro |103/88| 92/87 89/87 85/67 74/13 74/4 75/15 | 91/15 | 131/8 |152/1|137/37|106/73 | 1209/495 
Quiroga 119/88| 103/94| 99/103 94/57 86/10 89/4 88/8 | 103/13 | 148/3 |161/2| 169/15 |122/73 | 1 381/470 
San Jerónimo |109/89| 94/88 90/95 84/62 76/16 78/3 78/6 | 97/20 | 142/0|156/0| 157/7 |119/55 | 1280/441 
Erongarícuaro | 91/54 | 83/59 80/53 79/76 68/17 69/4 70/1| 80/22 | 115/0 |124/0| 123/9 |101/67|1083/372 
Cucuchucho |107/86| 95/86 91/84 84/62 75/13 74/4 73/9 | 94/16 | 137/0 |152/0|146/16|124/63 | 1252/439 
Granada 104/87| 90/84 86/90 80/63 72/18 74/3 73/7 | 94/21 | 139/0|154/0| 152/9 |118/48 | 1236/430 
Tzintzuntzan |114/93| 100/92| 95/100 90/62 82/13 84/4 84/8 | 100/18 | 145/3 |159/0|164/15|121/63 | 1 338/471 
Ojo de Agua |119/90| 104/95| 100/103 94/57 86/10 89/4 88/7 | 103/14 | 148/3 |161/2|170/13|123/72 | 1385/470 
Máximo 119/93| 104/95| 100/103 94/77 86/18 89/5 88/16 | 103/22 | 148/8 | 161/4|170/37|128/83 | 1 385/503 
Mínimo 91/46 | 83/54 80/41 79/43 68/8 69/3 70/6 | 76/13 | 106/0|119/0| 96/7 | 96/48 | 1083/308 
Media 109/81| 96/83 92/83 87/63 77/13 79/4 79/10 | 93/17 | 135/2|149/1|116/17|116/66 | 1261/440 
Desv. est. 8/17 7/15 6/21 6/10 6/3 7/1 6/3 9/3 14/3 | 15/1 | 18/10 | 11/10 97/60 


Nota: los valores en negritas correponden a la evaporación real. 


Cuadro 3. Repartición mensual de la evaporación del lago y la vegetación hidrófita. 


2006 2007 


Julio | Agosto|Septiembre | Octubre | Noviembre | Diciembre | Enero | Febrero | Marzo | Abril | Mayo | Junio Total 


Lago 99 85 81 74 66 68 67 Sil 136 151 145 116 110 
Lirio 66 66 64 66 64 66 66 59 66 64 66 64 777 
Thyía y tule 88 89 86 89 86 89 88 80 88 86 89 89 1 047 
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Elevación (m) 


Exageración vertical = 10x 
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Figura 2. Ubicación en altitud de los evaporímetros tipo ETgage y las estaciones EMAs. Los dos primeros perfiles 
corresponden a la dirección norte-sur (2-2” y 4-4) y los otros dos a la dirección este-oeste (1-1” y 3-37). 
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Figura 3. Variables climatológicas durante julio 2006-2007: temperatura, T; humedad relativa, Hr; 
déficit de presión de vapor, DPV; velocidad del viento, U,; y la radiación global y neta, R, y R,, 
registradas en las estaciones de Tupátaro (7), Pichátaro (O), Yotatiro (y), Chupícuaro (4), 

La Pacanda (m) y Pátzcuaro (0). 
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Figura 4. a) Comparación de la precipitación histórica acumulada del periodo 1969-1991 


en la estación Pátzcuaro, y la del periodo julio 2006-junio 2007 por subcuencas, incluyendo la del lago. 


Las cajas indican la media y mediana, y los percentiles del 5 y 95% de la lluvia histórica 1969-1991; 


b) distribución espacial de la precipitación del mes de junio 2007. 
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Figura 5. a) Comparación de la evaporación mensual estimada en la estación de Yotatiro, 
subcuenca de Ajuno versus Penman-Montheith-FAO, Turc, Turc (Abtew, 1996), Priestley-Taylor e Irmak; 
b) evaporación estimada para el lago de Pátzcuaro con los métodos de Priestley-Taylor, deBruin-Keijman, 
Penman (1948) y Makkink. Las cajas indican la media y mediana, y los percentiles del 5 y 95%. 
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Figura 6. a) Distribución espacial de la evapotranspiración (mm) de referencia calculada con el método 
Penman-Monteith-FAO para el mes de junio de 2007; b) distribución espacial de la 
evapotranspiración real, aplicando la corrección de Zeng (Sharif y Miller, 2006). 
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Figura 7. a) Comparación de la evaporación, estación Bowen versus método de Pennman (1948) y el ETgage; 
b) comparación de la evaporación ETgage versus Penman (1948); c) plano de exclusión de errores 
de la relación de $ versus Ae, en función de R,-G y distribución de los gradientes de presión de vapor (Ae) 
y de temperatura (AT) del aire registrada en la estación Bowen situada en la isla de la Pacanda. 
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vegetación hidrófita. El análisis del periodo 
julio 2006-junio 2007 refleja el estado actual 
de la repartición de la evaporación. En este 
sentido, la repartición de la evaporación a largo 
plazo puede cambiar en forma significativa 
debido a las alteraciones que está sufriendo 
la cuenca por el cambio acelerado del uso 
del suelo. El estudio enfatiza la metodología 
utilizada en este trabajo, las observaciones 
directas con el método BREB y el ETgage. 

Los resultados exponen que los métodos 
deben adaptarse a condiciones particulares de 
cada sitio para reducir las diferencias por los 
límites de aplicación, independientemente 
de las condiciones locales de las zonas. Por 
otra parte, no obstante que el año de estudio 
corresponde a uno relativamente seco, en 
comparación con la media histórica de la 
precipitación, existe un excedente para la 
generación de escurrimientos superficiales 


y pone 
en evidencia la importancia del lago en la 


e infiltración / drenaje profundo, 


regulación del flujo evaporativo que demanda 
la atmósfera dentro de la cuenca y el papel que 
juega en la regulación del clima dentro de la 
misma. 

La evaporación conjunta de la transpiración 
y la asociada con las pérdidas por intercepción 
oscila entre el 81 y el 88% de la precipitación, 
esto sin considerar la evaporación adicional y 
las pérdidas correspondientes a la vegetación 
herbácea y lechos vegetales, que en conjunto no 
son mayores al 7%, es decir, el agua disponible 
para los escurrimientos superficiales y la 
infiltración percolación es del 12-19% de la 
precipitación anual. 
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Abstract 


GONZÁLEZ-SOSA, E., MASTACHI-LOZA, C.A., RIVERA-VÁZQUEZ, J.B., GUTIÉ- 
RREZ-LÓPEZ, A., LAFRAGUA, J. £é GUEVARA-ESCOBAR, A. Evaporation in the 
Lake Patzcuaro watershed, Mexico. Water Technology and Sciences, formerly Hydraulic 
engineering in Mexico (in Spanish). Vol. L, no. 3, July-September, 2010, pp. 51-69. 


Annual evaporation in the Lake Patzcuaro watershed was  partitioned as 
follows: evapotranspiration of terrestrial vegetation, interception losses, lake evaporation 
and evapotranspiration of aquatic plants. Terrestrial evapotranspiration was estimated 
from the records of six stations and seven atmometers (ETgage). Lake evaporation was 
measured with the BREB method. The interception loss was calculated with lineal relations 
to different classes of vegetation, and evapotranspiration of aquatic vegetation with 
local measurements by open path chamber. The evapotranspiration of vegetation was 308 
mm <ETR <503 mm (+60), lake evaporation 1179 mm, aquatic evapotranspiration (ET H) 
1 824 mm, and interception losses amounted to 102 mm < El < 227 mm (+ 51). Water 
balance by subwatershed showed a deficit of 12 mm in Ajuno and over 65% in the lake. In 
the rest ofthe subwatershed the water balance was positive for surface runoffand infiltration 
or deep percolation: 18 mm <(P-TE <156 mm (+ 53). The results show the importance of 
vegetation and the lake in the distribution of evaporation. 


Keywords: evaporation, water balance, lake, BREB. 
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